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Abstract. We discuss recent Cu if -edge resonant inelastic scattering (RIXS) results 
for CUB2O4, a model system consisting of electronically-isolated CUO4 plaquettes, 
which exhibits anomalies in the incident photon energy dependence that are not 
captured by electron-only models. We show that these anomalous features can be 
qualitatively described by extending the electron-only considerations to include the 
lattice degrees of freedom. Our findings have important general implications for the 
interpretation of K-edge RIXS results for electronically higher-dimensional transition 
metal oxides. 
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1. Introduction 

Resonant inelastic X-ray scattering can be regarded as a resonant Raman process which 
results in electronically excited final states reached through intermediate core-hole states 
of X-ray edges. This species-specific, momentum-tunable probe is enabling the study 
of Mott-Hubbard physics [U El |3j H], Kondo phenomena in rare-earth intermetallics 
[5], d—td excitations of transition metal species in complex oxides [6j [7] and, more 
recently, orbital [8] and magnetic [9J [10] excitations. While this promising technique 
has led to new insight into the excitation spectra of complex materials, its status as 
an experimental probe is relatively immature in comparison to the fields of optical or 
photoemission spectroscopy. 

Here we describe the consequences of common considerations [3, [12] of the RIXS 
cross section of a single C11O4 plaquette, resonantly excited through the Cu K (Is— )4p) 
X-ray absorption edge. Motivated by recent experimental work for the relatively simple 
model system CuB 2 04, we then point out aspects of experimental data pj, HHJ [T6] 
which cannot be explained by electron-only considerations [12]. The discrepancies 
can be directly addressed through consideration of the coupling of lattice degrees of 
freedom through the standard Franck-Condon formalism. We conclude that similar 
considerations are necessary in the analysis of X-edge RIXS results of more complex 
transition metal oxides. 

2. Electron-only Model 

With only one hole, the ground state of the single CUO4 plaquette transforms as x 2 — y 2 , 
with a completely filled Is electronic core shell, a completely empty Cu Ap shell, and one 
hole shared among the Cu 3d and O 2p orbitals[llJ. A photon with sufficient energy (~ 
8995 eV) can be absorbed, inducing a dipole- allowed optical transition, which results 
in a Ap electron and a Is core hole. In the simplest treatment of Cu K-edge RIXS, the 
influence of the core hole is to locally repel valence holes, thereby disturbing the Cu 
3d and O 2p valence subsystem, while the extended Ap electron is considered to have 
a relatively limited influence on the dynamics [HJ [71 [321 CEl IHj. The RIXS process is 
completed when the Ap electron radiatively recombines with the Is core hole, emitting 
a photon which leaves behind an amount of energy reflective of a core-hole-induced 
excitation created in the valence system. Because the core hole potential is localized, 
the resulting spectroscopic process results in electronic transitions among valence states 
of the same symmetry, so we may limit our attention to states with x 2 — y 2 character, 
as is commonly done in considerations of core-level photoemission spectroscopy [T3j . 

Following van Veenendaal, Eskes, and Sawatzky [13], we therefore need only 
consider the x 2 — y 2 state \d 9 ) with hole density on the Cu and on the O 2p ligand 
state of the same symmetry, which we call \d 10 L) (sometimes called \p x 2_ y 2) [T3]). 
The wavefunctions for these states are shown in Fig. la. These basis states are 
energetically separated by the bare charge-transfer energy A, but are further split under 
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Figure 1. (a) Hilbert subspace of a single CUO4 plaquette. Basis states (d 9 and 
d 10 L) are shown as black arrows. The effect of introducing hybridization t vd and 
rediagonalization is effectively a rotation in this picture, giving the red arrows, which 
show the ground (g) and molecular orbital (MO) states. Further introduction of the 
core hole potential gives the purple arrows, depicting the well-screened (WS) and 
poorly screened (PS) states. The matrix elements for various spectroscopic processes 
are represented by dotted lines. Hole wave functions are shown around the perimeter, 
and the angles and <p used in the text are defined, (b) Energy level diagram for this 
system, showing the effect of including t p d and the core hole potential. Colors used 
correspond to (a), (c) RIXS profile of the single plaquette, after [7]. The state energies 
are: E g 0, E MO 6 eV, E ws 8990 eV, E PS = 8998 eV. 

the hybridization t p d, and the new ground state is formed as a bonding combination 
of the two basis states: \g) = cos 9\d 10 L) + sin#|d 9 ). The antibonding state \MO) = 
— sin 9\d 10 L) +cos 9\d 9 ) is generally referred to as the molecular orbital (MO) excitation. 
Transitions from the ground to the MO state are RlXS-active at the Cu K edge because 
they do not involve changes in symmetry of the electronic wavefunction [121 03] • This 
is very different from the case of soft X-ray RIXS at the Cu L 3 edge j6j EJ HUJ [T7] . 
where symmetry-changing d^d transitions are favored over these symmetry-preserving 
"charge-transfer" excitations. 

During the Cu K-edge RIXS process, a deeply bound Is electron makes a dipole- 
allowed transition to a spatially extended, high-energy Ap state, leaving behind a Is 
core hole. This strong perturbation affects the structure of the valence electronic states, 
giving rise to two new states, shown as purple arrows in Fig. la. The lower-energy "well- 
screened" state \ WS) = cos cj)\d 10 L) +sin0|rf 9 ) has a large electron density screening the 
core hole potential, while the "poorly-screened" state \PS) = — sin <j)\d 10 L) + cos0|d 9 ) 
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has a higher hole density. In Cu L 3 X-ray photoemission spectroscopy (XPS) [13J, 
the projection of the ground state onto each of these states determines the spectral 
profile (e.g. Jxps oc Kl^S^t/)! 2 ). However, in a resonant elastic or inelastic scattering 
experiment (REXS or RIXS), the signal arises from double projections, first from the 
ground state (g) onto highly excited intermediate states (WS or PS), and then back to 
the lower-energy excited state (MO). In our case, this double projection has a particular 
property: (MO\WS)(WS\g) = -(MO\PS)(PS\g), independent of t pd or core hole 
interaction strength, as can be seen graphically in Fig la. This fact implies that there 
is equal strength of RIXS scattering through each of the two intermediate states into 
the MO final state, leading to a double-peak pattern (Fig. lc). A streak of intensity 
between these two resonances occurs due to constructive interference for intermediate 
photon energies, and is a consequence of the resonant Raman process for two non- 
degenerate intermediate states with equal and opposite matrix elements [7j. 

3. MO excitation in experiments 

To date, the MO excitation has been observed in cuprates of various electronic 
dimensionalities: quasi-zero dimensional [7], edge-shared one-dimensional [T2l [16], 
corner-shared one-dimensional [18] , two-dimensional planar (TJ [2j CEH1 |3j [4] , and three 
dimensional [20]. The electronic structure is far more complicated in the higher- 
dimensional systems than in the molecular limit, but a localized MO excitation is 
nonetheless observed [T5] . As discussed above, it appears that this excitation resonates 
strongly at both the WS and PS states of the low-dimensional systems H2]. In 
general, an approximate correspondence can be made with the WS and PS states 
discussed above, and it is through a higher-energy PS-like state that the MO excitation 
resonates most strongly, due to the localized nature of the dominant d 9 character. A 
WS-like intermediate state can also be identified, but because this state involves the 
displacement of hole density away from the Cu site and onto the neighboring O, it 
couples most directly to extended excitations of both the charge-transfer [2j [TJ [191 El H] 
and magnetic [9j [10] types. 

Figure 2a shows Cu K-edge RIXS data for electronically quasi-zero-dimensional 
CuB 2 04. The data were collected at SPring-8 BL11XU with instrumental energy 
resolution 380 meV FWHM. The strong Raman-dispersing peak at ^ 7 eV transfer is the 
MO excitation [7J . The low-energy d—td and higher-energy 4p excitations are identified 
and discussed in detail elsewhere [7J. Fig 2b shows an anomalous 20% sharpening of 
the excitation at an incident energy 9002 eV. Figure 2c shows that there is a clear 
step-like incident-energy dependence to the MO feature's energy, which mimics the 
behavior found in the edge-shared cuprate LiCu 2 02 [32], also shown in Fig. 2c. These 
observations can not be understood in the electron-only treatment discussed above, and 
were not reproduced by the purely electronic cluster calculations of Hancock et al [7] or 
Vernay et al. [12]. A third puzzling aspect of the data is the broad Gaussian lineshape, 
which is unexpected for a long-lived electronic final state that is weakly coupled to its 
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environment. We interpret these salient, but enigmatic properties to imply that purely 
electronic considerations are insufficient to fully understand the data. 
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Figure 2. (a) Line scans for CUB2O4 [7] showing the incident-energy dependence 
of the MO excitation, (b) Incident-energy dependence of the MO Gaussian width. 
The line shows the same for the model discussed in the text, multiplied by a factor 
of 2.5. (c) Incident-energy dependence of the MO peak position of CUB2O4 (incident 
polarization e^|| (001)) [7 (solid circles) and LiCu202 [E]. The line shows the position 
from the model calculation (Fig. 3b). (d) The incident-energy dependence of the 
Gaussian position for La2Cu04 (from [1]). 



4. Effect of lattice degrees of freedom 

We now make the simplest possible extension of the above model to include the lattice 
degrees of freedom of the system, and show that these three puzzling aspects of the data 
can be succinctly addressed. The procedure follows the Franck-Condon (FC) treatment 
for a molecular system [2T1 [T71 [22] . wherein we consider as before transitions among 
the electronic degrees of freedom (e), but with additional coupling to the vibrational 
continuum (v) within the sudden approximation at each stage of the two-step RIXS 
process. This process is shown schematically in Fig. 3a. 

The RIXS cross section for this model is then given by the Kramers-Heisenberg 
formula: 

da 00 

— oc \M\ 2 S(h(u in - v sc ) - (E M o + niMohuo - E g )). (1) 

d ^ 1 m=0 
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with the amplitude 

M = y y [(MO\i)(i\g)] e [(m M o|^)(n,|0,)], 

Here ujq is the natural frequency of the harmonic potential, T is the core-hole decay rate, 
Emq, Ei, and E g are the (electronic) energies of the MO, intermediate (WS or PS), and 
ground states, respectively, and fiv in (hv sc ) is the energy of the incident (scattered) 
photon. 

The result for the purely electronic model in Fig. lc is recovered if we assume 
negligible lattice coupling and molecular potential deformations. As before, we would 
have a strong MO excitation, resonating equally at the WS and PS states, whose position 
is given by the bare electronic transition energy E M o — E g , and whose width is narrow, 
determined by the intrinsic lifetime of the final electronic state. The new elements are 
the Franck-Condon factors, which represent the overlap of nuclear wavefunctions from 
the ground and excited state potential energy surfaces. These overlaps can be calculated 
analytically [22] if we assume a purely translational shift of a one-dimensional harmonic 
potential. The FC factors are then calculated as simple overlaps of harmonic wave 
functions (Fig. 3a), e.g.: 

(a \n-N / \ \ n+m-2N 

[{m a \n b )\ v = e W - - _ — L N { %a ) (3) 

where N = min(n,m), is the associated Laguerre polynomial, and 7^5 = 
\/^h( Xb ~ Xa ) * s ^ e electronic transition-induced displacement of the harmonic potential 
in units of the ground state positional uncertainty. 

The parameters of the model which express the influence of phonon coupling are 
r ) g ^i ) and 7^mo, where i G {WS, PS}, and the phonon frequency uj . The latter sets the 
energy scale for coarse changes in both directions of the energy transfer/incident energy 
plane. In this plane, the core-hole broadening (r~l eV) provides a large electronic 
energy uncertainty in the incident energy direction. Therefore, the summation over 
the intermediate state phonon index n has highly overlapping, interfering contributions, 
with relative phases which are energy dependent close to the resonance. The result of 
this summation also depends on the final state phonon index m, and therefore weakly 
affects the lineshape in the energy-transfer direction. The model therefore supports 
incident-energy dependent lineshape properties, with even one electronic transition, in 
contrast to purely electronic models [12] . 

In more detail, the lineshape is dependent on the relationship between 7^ and 
Ii^mo- If Ig^i+li^f = 7^^/ = 0, there is no net displacement of the potential energy 
surface and no effect in the energy transfer direction: we again recover Fig. lc, with an 
expected narrow Lorentzian. For a fixed 7#^, moving 7^/ in either direction away from 
the symmetry point 7^/=0, the lineshape blueshifts and quickly becomes Gaussian for 
even the modest value of |7^/| ~ 3. Figs. 3c and 3d respectively show the position 
and width for a range of values of 7^/ with 7^ = +3. In addition to the blueshift of 
order fiuj^^f and a broadening of order /kj |7^/|, step-like modulations appear in the 
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Figure 3. (a) Schematic of the Franck- Condon process. The electronic transitions 
are accompanied by a manifold of vibrational transitions, so that overlaps of the 
nuclear wavefunctions between ground, intermediate and final states need to be 
considered. The overlaps strongly depend on the relative displacement of these three 
states (coordinates x 9 , xws, xps, an d #mo)> which enter the FC amplitudes as 
ja^b — \J ^j%{ x b — x a ). The one-dimensional spatial ordering of these four states 
is implied by the fit to the model shown in Fig. 2c. (b) Plot of Eq. (1) for the FC 
factors given in Eq. (2) with j 9 ^ W s -4, Jws^mo 8, Ig^PS 6, and Jps^mo 
-2, with hujQ = 75 meV. The vertical line at 5.55 eV shows the bare electronic excitation 
energy. Incident energy-dependent (c) position and (d) width for fixed Jg^i — +3 and 
various values of 7^/. 

energy position and width when the incident photon energy passes through resonance. 
The direction of the step in incident energy position depends on whether < 7^/7^ 
(region A) or > 7^/7^ (region B). These two limits are sketched in Fig 3a as a set 
of red or blue arrows, respectively. 

The comparison in Fig. 2b of these trends in the incident-energy-dependent position 
in the data for CuB 2 04 and the calculation suggests that the PS state is in region A. This 
assignment is consistent with the fact that there is an abrupt stiffening and narrowing 
of the feature as his in is passes through the PS resonance around 8998 eV. The situation 
for the WS state is less clear, due to the lack of narrowing at lower incident energies, 
as well as likely interference from the approaching Kfa^ fluorescence line, which is a 
relatively strong 3d^ls radiative transition at fixed scattered photon energy 8976 eV 
(Fig. 2a). Despite ambiguities regarding the description of the WS state, we tentatively 
put the WS state in region B based on the step in the MO position near 8994 eV. Fig. 
3b shows the model results with E PS E ws =5-55 eV and ("jg^ws, Iws^mo, Ig^ps, 
Ips^mo) = (-4, 8, 6, -2), corresponding to 7^/=4. The vibrational energy was chosen 
to be huoo = 75 meV, which is a typical value for bond-stretching phonons in cuprates 
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[23] . The values of the four parameters 7 imply a spatial ordering to the equilibrium 
positions shown in the schematic of Fig. 3a. 

Having established that the phonon-coupled model in Eqns. ((l)-(3)) is capable of 
describing the data at a qualitative level, we point out limitations that may be important 
to a quantitative application. First, the experimental value of the MO width (<~ 760 
meV) is about 2.5 times larger than the model value. The quantitative comparison has 
several potential shortcomings. One is the large uncertainty in the intrinsic width of 
the electronic excitation, which could strongly effect a quantitative analysis based on 
the widths of the MO feature. If we assume a larger intrinsic width, we can partially 
remedy the discrepancy, but still underestimate the narrowing near fiv in = 9000 eV. 
Another point is that the curvature of the internuclear potential energy is not likely 
to be the same at each stage of the process, as we have assumed. Such effects would 
manifest themselves as different values of u)q for the g, MO, WS, and PS states. In an 
extreme case, the potential energy surfaces may not even support bound vibrational 
states, but rather have an unstable, dissociative nature. Additional complications arise 
from the fact that CuB 2 04 contains two crystallographically distinct plaquette types, 
which may inhomogeneously broaden the lineshape, leading to an overestimate of the 
importance of phonons based on the width alone. Furthermore, there exist scattering 
channels which are not included in the above considerations, such as those which give 
rise to the d^d and transitions described in [7j. It is therefore possible that the 
peak is composed of several electronic transitions, further obfuscating an interpretation 
of its width. While a detailed quantitative application of the FC model with quantum 
chemistry methods is beyond the scope of the present work, we have shown that this 
model has the desired properties to describe the Gaussian lineshape with concurrent 
resonance-induced narrowing and step in energy transfer for the PS state. 

5. Implications, conclusions, and outlook 

We now discuss some general implications of our observations. Kim et al. [15] performed 
a systematic study of the molecular orbital excitation energy for a wide range of cuprates, 
and obtained the surprising result that this energy scales with the copper-oxygen bond 
length rfcu-o as d(5u-o with a ~ 8±2. This behavior is very different from that 
expected based on Harrison's rules [24] and band structure calculations [25], for which 
the exponent is in the range a^3.5-4. We identify the feature observed in RIXS as 
the maximum of a Franck-Condon band (at energy E M o + ^ohg^fl ~ Eg)i an d not 
as the bare electronic transition (zero-phonon line at E M o — E g ), which would reside 
at much lower energy (e.g., the vertical line in Fig. 3b). Consequently, in addition to 
a purely electronic contribution with small exponent, we expect a significant additive 
contribution that arises from the bond-length (dcu-o) dependence of the molecular 
potential shifts (7^/ and 7^). In our interpretation, this dependence is determined by 
the way the relevant phonon frequencies and potential shifts depend on the bond length. 
Our result suggests that the parameter combination fkj |7^/| behaves approximately 
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Reference [7j assigns the weak ^2 eV of CuB 2 04 (reproduced in Fig. 2a) to a d^d 
excitation and showed that this transition is a robust, one-plaquette effect, and not the 
result of local symmetry-lowering distortions. This transition is anomalous, because it 
necessarily involves a change in symmetry, an effect unexpected at the Cu K edge. In 
the context of phonon coupling, we can address its appearance by considering the effect 
of the Ap electron on the molecular potential surface. If the photoelectron reaches, say, 
a Ap x state, anisotropic screening could force the magnitude of molecular displacement 
of the O placed along the ±x axis to be different from the displacement of the O placed 
along the ±y axis. This distortion of the plaquette in the B\ Q channel in the intermediate 
state could open a pathway to a d^ z 2_ r 2 final state, for example. In Ref. [7], it was found 
that the d^d excitation resonates more strongly at the WS and PS states with Ap Xjy 
polarization than with Ap z polarization, thus supporting this scenario for d—±d genesis. 
Future high-resolution data for CuB 2 04 or other low-dimensional cuprate systems may 
help give better insight into the nature of the d^d transition and the plausibility of this 
scenario. 

Finally, we discuss the possible extension of the lattice-coupled RIXS cross section 
to planar cuprates. Abbamonte et al [lj reported a 6 eV excitation at the K edge 
of Cu in La 2 Cu04, which was described by a Gaussian line shape dressed by resonant 
denominators representing the incident and scattered photon poles. Later work in a 
different scattering geometry [HI El H] resolved a fine structure of excitations, so the 
original work in Ref. [lj should be regarded a measure of the overall inelastic spectral 
weight, containing both the MO excitation as well as charge-transfer excitations. As 
seen from Fig. 2, the incident-energy-dependent position of this spectral weight exhibits 
qualitative similarities with the MO excitation results for CuB 2 04 and LiCu 2 2 as well 
as our model calculation. The magnitude of the shift is much larger for La 2 Cu04, 
perhaps reflecting the enhanced range of influence of the core hole potential in generating 
phonons, or additional effects due to the complex electronic structure of the Cu0 2 plane. 
More recent experimentation [U [2j [191 HI EE] revealed that a sharp, relatively weak 
2.2 eV charge-transfer feature resonates at the WS state of La 2 Cu04 and appears to 
have a temperature-dependent lineshape [26]. The latter observation may be addressed 
already within a simple extension of our model: we have considered transitions only from 
the electronic- vibrational ground state, but transitions from the ensemble of thermally 
activated vibrational states would inevitably change the line profile, suggesting a route 
to understanding the observed temperature dependence at the relatively high energy of 
- 2.2 eV. 

In conclusion, we have shown that the electron-only models are incapable of 
describing even salient experimental results at the transition metal K-edge. On the 
other hand, we demonstrate for the molecular-orbital excitation of the simple model 
system CuB 2 04 that considerations of the coupling of the lattice to electronic degrees 
of freedom leads to a natural interpretation of the observed line shapes and (incident) 
photon-energy dependence. Nonadiabatic effects were first observed in C K-edge RIXS 
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of ethylene [17] . and our results show that these effects carry to the 3d transition metal 
edges as well. Phonon coupling is likely present in existing RIXS data for cuprates 
111 E | and other transition metal oxides [HI [27], as we have shown for LiCu 2 02 
[T5] and La 2 Cu04 pQ. We have also identified an important property of the model 
that the lineshape away from resonance is most heavily influenced by the electron- 
phonon coupling between the ground and final state-it is only near resonance that 
the intermediate state effects play a role. Finally, we emphasize that the process 
discussed here is only one possible channel for RIXS scattering, where the photoelectron 
recombines, and the charge excitations occur only via the core hole interaction. In the 
case of d^d excitations [6], orbiton scattering [8], and magnetic [9j[l0] excitations, the 
applicability of the present considerations is not as clear. We expect future experiments 
with higher energy resolution, coupled with more refined theoretical considerations, 
to directly address the role of the electron-phonon interaction in determining RIXS 
lineshapes. 
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support of Eleiza Braun. This work was supported by the DOE under Contract No. 
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